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Abstract: We explore as clearly as possible the features of neutrino oscillation 
which are relevant for measurements of the CP violating Kobayashi-Maskawa phase 
6 and the sign of Amf^. We focus on the so called low-energy option and discuss 
principles for optimizing experimental parameters to measure these two quantities 
simultaneously. Toward the goal, we first formulate a method for obtaining a bird-eye 
view of the phenomenon of neutrino oscillation by introducing a new powerful tool 
called the "CP trajectory diagram in bi-probability space". It allows us to represent 
pictorially the three effects separately in a single diagram; effect from genuine CP 
violation due to the sin 5 term, effect from the CP conserving cos 5 term, and the 
fake CP violating effect due to earth matter. By using the CP trajectory diagram we 
observe that there is a two-fold ambiguity in the determination of 6 which is related 
with the sign of Amfg. We then address the question of what are the promising 
options for conceptual design of experiments at low energies which looks for CP 
violation and at the same time would resolve the two-fold ambiguity. We point out 
that a version with distance of about 700 km, CERN to Gran Sasso and/or Fermilab 
to Soudan-2 site, with a megaton class water Cherenkov detector gives an optimal 
design which allows simultaneous determination of 6 and the sign of Am^g in situ. We 
also point out that there is a possibility that the similar in situ measurement of both 
quantities can be done at the Phase II of JHF experiment with much shorter baseline, 
under the assumption of nature's kind setting of 6 to the region of sin 6 ■ Amf^ < 0. 
A technique of running at high (~ 1 GeV) and low (~ 0.5 GeV) beam energies is 
proposed as a method for better identification of 6. 
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1. Introduction 

Discovery of neutrino oscillation in atmospheric neutrino observation in the Super- 
Kamiokande (SK) experiment |l| opened up a new window to physics beyond the 
standard model of particle physics. Moreover, robust discrepancy of the measured 
flux of solar neutrinos to the calculated one 0] presents another indication for neu- 
trino masses and the lepton flavor mixing. In fact, the flrst result from the SNO 
group combined with the results by the Super-Kamiokande solar neutrino obser- 
vation 0] strongly indicates that solar neutrinos also do oscillate. The existence 
of phenomenon of neutrino oscillation is further strengthened by the result of the 
K2K experiment, the flrst long-baseline experiment with artiflcial neutrino beam, in 
particular by their latest result 
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The determination of the complete structure of the lepton flavor mixing matrix, 
the Maki-Nakagawa-Sakata (MNS) matrix is one of the most challenging task 
in particle physics. While we began to grasp the values of the leptonic mixing 
parameters there remain three quantities which are poorly known or not constrained 
at all at this moment. They are ^13, the sign of Am^g, and 6, the leptonic Kobayashi- 
Maskawa angle 0, ||. (See e.g. Ref. || for a summary of situations of yet to be 
determined parameters in the three-flavor mixing scheme of neutrinos.) Regarding to 
6'i3, we only know its upper limit, sin^26'i3 < 0.1, from the reactor experiments [p!0|] . 
The bound is to be improved, or the value of 613 itself could be determined by the 
next generation long baseline experiments [0, [1^, [1^. About the sign of Am^g, while 
there is a strong indication that the inverted mass hierarchy (or negative Am^g) is 



disfavored by the observed neutrino events coming from supernova SN1987A [|M|, no 
hint is available from laboratory experiments. With respect to CP phase, 6, we do 
not have any experimental clue at all. 

In this paper, we intend to explore as clearly as possible the features of neu- 
trino oscillation which are relevant for experimental measurements of the latter two 
quantities, 6 and the sign of Am^g. We focus on the so called low-energy option 
and discuss principles for optimizing experimental parameters to measure these two 
quantities simultaneously. Toward the goal, we flrst formulate a method for obtaining 
a bird-eye view of the phenomenon of neutrino oscillation at low neutrino energies, 
typically, 0.5-2.0 GeV. We do this by introducing a new powerful tool called the "CP 
trajectory diagram in bi-probability space", which allows us to represent pictorially 
the effects of genuine CP violating phase and the earth matter separately in a 
single diagram, as explained in detail below. We then address the question of what 
are the promising options for conceptual design of experiments which look for CP 
violation among thinkable varying possibilities at low energies. 

We have discussed in previous communications [|TB], [I^ the possibility of use of 
an intense low energy neutrino beam, which is nowadays referred to as a "super- 
beam" a la Richter |TB[, to measure CP violation. It is based on the underlying 
cancellation mechanism of the leading-order matter effect, the vacuum mimicking 



mechanism 19 , which would allow us to measure CP violation in a clean vacuum- 



effect-dominated environment for neutrino oscillations. It gives us a great merit of 



avoiding the notorious problem of matter effect contamination pO|, ^ |2^ . 

As a concrete example which realizes our basic idea we designed and examined 
an experiment which uses a neutrino superbeam of energy E ~ 100 MeV and a 1 
megaton water Cherenkov detector |]T6|, [l^. We have ended up with a rather short 
baseline, L ~ 30 — 40 km, as optimal distance in such an experiment. While it served 
well for our purpose of illuminating the basic idea, several experimental problems 
have been raised when it was taken at face value as an experimental proposal. (See, 
however, for a detailed feasibility study of the very similar idea, E = 250 MeV, 
L ~ 100 km, and a 40 kton water Cherenkov or liquid scintillator detectors.) 

If we simply scale up the energy range to the more realistic one, E 1 GeV, 



keeping E/L flxed as in Ref. |]I6|, [T^, the optimal distance may become naively L ~ 
300-400 km. It is very similar to the experimental set up recently discussed by the 
JHF neutrino experimentalists. In particular, by the JHF neutrino working group. 
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various options for neutrino beam, wide band (WB) beam, narrow band (NB) beam, 
and off-axis (OA) beam with neutrino energies of 500 MeV-a few GeV have been 
extensively studied p9|. Moreover, the efficiency of removing vr^ contamination has 



been improved tremendously since the old Letter of Intent (LOI) JlTI by implement- 
ing severe cut by imposing a second ring, a new technique originally developed by 
the Super-Kamiokande group [^]. See their new version of LOI in 

In this paper we try to make a step forward along the line of thought toward 
measuring CP violation with use of low energy conventional superbeam by describing 
a general strategy of optimizing beam energy and/or baseline distance. We will show 
that the matter effect are comfortably large even in these mediumly long baseline 
(~ 300 km) experiments so that there is a possibility of simultaneous measurement 
of 6, the leptonic Kobayashi-Maskawa angle, and the sign of Aml^f.^ in situ in a single 
experiment. 

An alternative strategy based on intense neutrino beam from a muon storage ring 



called neutrino factory has been extensively discussed in the literature p5|, |2g, ^ 
Our strategy which utilizes low energy conventional superbeam differs in many ways 
from neutrino factory, e.g., on beam energy, detection principle, and most crucially, 
on whether the value of ^13 must be known in advance (conventional beam), or can 
be measured simultaneously with 6 (neutrino factory). While many discussions are 
going on upon ignition by Ref. about which strategy is more superior [P^], 

we strongly feel that both strategies must be fully developed both in physics and 
beam technology aspects before attempting any real comparisons between the two 
strategies. 

The motivation for our consideration of neutrino experiments with beam energies 
of ~ 1 GeV in this paper is partly theoretical and partly experimental. Experimen- 
tally, the energy region comes out as a natural compromise of the two requests that 
CP violation is large and neutrino beam is intense enough. Furthermore, it has a 
great merit of being able to utilize the results of the recent developments we just men- 
tioned above. On the theoretical side, we discuss below by using the "CP trajectory 
diagram" a principle of optimizing beam energy and/or path length to maximizes 
the detection probability of CP violation. We will see that such discussion naturally 
leads to several options which utilize the same energy region, ~ 1 GeV. Notably, 
we will uncover the possibility in which an in situ simultaneous measurement of CP 
violating phase 6 and the sign of Am^g will be possible for (i) whole region of 6 
{L ~ 700 km), and (ii) half a region of 6 which fulfills the condition sin 5 ■ Amfg < 
(L ~ 300 km). 

In Sec. 2 we point out that an approximate two-fold degeneracy exists in vacuum 
neutrino oscillation probability and show that the degeneracy is partly lifted by the 
matter effect. We introduce in Sec. 3 the CP trajectory diagram on bi-probability 
space, and thoroughly analyze its properties. It will be shown that it is a powerful 
tool for illuminating the general structure of the oscillation probabilities of neutrinos 
and antineutrinos. We will point out that the two-fold degeneracy which exists 
in vacuum neutrino oscillation probability, after resolved by matter effect, leaves a 
remnant ambiguity in the determination of CP violating phase 6. We then discuss 
in Sec. 4 a principle of tuning beam energies to have maximal CP violation as well 
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as to help in resolving the two-fold ambiguity. In Sec. 5 we describe several possible 
ideas for resolving the two-fold ambiguity. Throughout these sections we will reveal 
that a new strategy toward simultaneous determination of the CP violating angle 6 
and the sign of Amfg naturally emerges from the discussions of problems mentioned 
above. In Sec. 6 we describe some concrete examples of experiments which utilize 
conventional neutrino superbeams. We estimate the number of events and include 
the background rate to obtain a rough idea for the accuracy of the measurement. In 
Sec. 7 we briefly discuss the CP trajectory diagram with experimental parameters 
for neutrino factory. In Sec. 8 we state our conclusions. 



2. Matter effect helps to resolve two-fold ambiguity in vacuum 
neutrino oscillation 



We start by describing a role played by matter effect to help resolve a two-fold 
ambiguity which would exist in a vacuum-effect dominated neutrino oscillation ex- 
periment to measure CP violating angle 6. It contrasts to the negative role played 
by the matter effect as a contamination in measurement of genuine CP violation, the 



widely recognized fact in the literature pO| , [2^ , |23[] . We hope that the discussion 
illuminates the necessity of complete understanding of the interplay between effects 
due to the Kobayashi-Maskawa phase and earth matter. 

We use throughout this paper the standard notation of the MNS matrix: 



U 



C12C13 



■S12C13 



~'5l2C23 ~ Ci2S23Si3e*'^ C12C23 — Si2S23'5l3e*'^ ■523Cl3 
S12S23 ~ Ci2C23Si3e*'' — C12S23 — Si2C23Si3e*'' C23C13 



:2.ii 



— ■■-J where Am^j 

is assumed to be smallest, relevant for oscillation solutions to the solar neutrino 



We define the neutrino mass-squared difference as Am'f- = nij 



problem. We note in passing that in the MSW mechanism, under the assumption 
that Am23 ~ Amfg is much larger than Am^g; latter can be made always positive 
as far as 9i2 is taken in its full range [0, 7r/2][^. 

The ambiguity arises due to the fact that Am^ scale implied by solar neutrino 
solution may be small such that 



Am^2-^ 
2E 



2.54 X 10 



1-2 



10-5 eV^ i (lOOO km) (l Gev) 



(2.2) 



is sufficiently small compared to unity for most of the possible experimental parame- 
ters. To first order in the parameter the neutrino oscillation probabilities of z/^ —>■ v^. 
and z/^ — > Ve in vacuum are given by 



sin^ 2^135^3 sin^ ' 



\ AE 



1 



' Am^2-^' 
2E 



sm 



2E 



4 



T4J,=in«(^jri„^(^j, (2.3) 

where — and + sign in front of the 4th term refer to the neutrino and the antineutrino 
channels, respectively, and Jr = £125120^3513023523 denotes the reduced Jarlskog fac- 
tor. One can readily observe that the oscillation probabilities, except for the second 
term, are invariant under simultaneous transformation 

6 ^ 71 — 6 (mod.27r), 

Amj, ^ -Ami,. (2-4) 

Under the transformation ( p.4|) , flipping sign of the cos 5 term by the first transfor- 
mation is canceled by the second, whereas the sin 5 term is manifestly invariant. It 
implies that the probability is approximately degenerate for two values of 6, unless 
one know a priori the sign of Ami,, hence there is a two-fold ambiguity in 
determination of 6. 

We note that the invariance holds in a very good approximation. It is because 
the second term is small in magnitude because of the suppression factors including 
^ Am^L ^j ^^^^ t\]_a,t sin^ 26'i3 is smaller than 0.1 in order to satisfy the CHOOZ 

result |T0[. Also the higher order terms of (— is sufficiently small for most of the 
experimental settings. It may be worthwhile, however, to note that the degeneracy 
is accidental and approximate in nature, and represents neither inherent nor exact 
properties of vacuum neutrino oscillations. 

Of course, the ambiguity does not exist if we know in advance the sign of Amf,. 
However, it is very unlikely that the question of the normal vs. inverted mass hi- 
erarchies of neutrinos will be answered in a convincing way in the near future. It 
is because (as it is believed) the determination of the sign of Amf, requires mea- 
surement of interference between the CP and the matter effects, which necessitates 
a sufficiently long baseline. It is not an easy experiment to carry out because it 
requires either intense neutrino beam or supermassive detectors, or plausibly both. 

We now point out that the two-fold degeneracy in vacuum oscillation probability 
is lifted by the matter effect represented by the index of refraction a{x) = y^GpNe^x) 
where Nf.{x) is the electron number density in the earth. When we include the matter 
effect there arise, in leading order of aL, the following additional terms Pmatt in the 
oscillation probability computed under the adiabatic approximation [2U|: 



r / N / M ^9^9 f 2Ea(x)\ 9 /Am?oL\ 

Pmatt±[^t^{^t,) M^e)] = ± COS 2^i3 slu' 2^135^3 I 1 slu' I j 

T sm^ 2^13 COS 2^135^3 sm I \ , (2.5) 

In ( |2.5| ), a{x) = \/2GFNe{x) where Gi? is the Fermi constant and Ne{x) denotes the 
electron number density at x in the earth, and +(— ) in the 1st term and — (+) in 
the 2nd term refer to the neutrino (antineutrino) channel. 
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It is easy to observe that the degeneracy is hfted because Pmatt is not invariant 
under flipping sign of Am^g, a well known fact. What is perhaps not so well known is 
that the lifting of the degeneracy does not completely resolve the two-fold ambiguity 
in determination of 6. We introduce in the next section a powerful tool called the 
"CP trajectory diagram" and demonstrate that a remnant ambiguity exists in the 
determination of CP violating phase. Since the ambiguity is related with the sign of 
Am^3, we are naturally invited to the problem of simultaneous determination of 6 
and Amfg, which we will pursue in Sec. 6. 

3. CP trajectory diagram in bi-probability space 

To illuminate global features of neutrino oscillations relevant for low energy experi- 
ments, we introduce the CP trajectory diagram in bi-probability space spanned by 
P(z/) and -P(z^) Unless otherwise stated we simply denote P{i^^ — > i^e) and 

P{'^fi i^e) as -P(z^) and P{i^) in this paper. We show in this section that the di- 
agram is a useful tool for our purpose because it can display pictorially the three 
effects separately in a single diagram; genuine effect of CP violating phase S coming 
from sin 5 term, CP conserving effect due to cos 5 term, and the matter effect. 

Suppose that we compute the oscillation probability -P(z^) and -P(z^) with a 
given set of oscillation and experimental parameters. Then, we draw a dot on two- 
dimensional plane spanned by -P(z^) and -P(z^). When 6 is varied we have a set of dots 
which forms a closed trajectory, closed because the probability must be a periodic 
function of 5, a phase variable. 

In Fig. 1 plotted is the contours of oscillation probabilities P(i^) and -P(z^) which 
are drawn by varying the CP violating phase, S, from to 27r. We note that in 
this work, while we refer to analytic expressions for the explanations, all the results 
shown in our plots as well as in tables were based on the computations obtained by 
numerically solving the three flavor neutrino evolution equation assuming constant 
matter density. The solid and the dashed lines are for positive and negative Amfg, 
respectively. The dotted and the dash-dotted curves are the cases in vacuum in 
which the matter effect is switched off in the corresponding cases of the solid and 
the dashed lines, respectively. In fact, the abscissa and the ordinate are not quite 
the oscillation probabilities but are the ones averaged over an appropriate energy 
distribution of (anti-) neutrinos. It is to avoid accidental zeros, and at the same 
time is also meant to mimic the average over the energy dependent flux times cross 
sections, the procedure we will actually execute in Sec. 6. The neutrino energy 
distribution is taken to be a Gaussian shape with central value of (a) 0.5 GeV, 
(b) 1.0 GeV, (c) 1.5 GeV, and (d) 2.0 GeV, respectively. The widths of Gaussian 
distribution is taken to be 20 % of the peak energies, e.g., 100 MeV for E = 500 
MeV in Fig. la. The baseline length is taken to be 295 km, JHF-Kamioka distance, 
and the values of other parameters are typical ones for the large angle MSW solution 
to the solar neutrino problem and sin^2^i3 = 0.05, as given in the caption of Fig.l. 
Notice that unless the large angle MSW solution is the case the measurement of CP 
violation would be very difficult. 
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Figure 1: CP trajectory in the bi-probability (given in %) plane for the baseline L = 295 
km. As indicated in the figures, the sohd and the dashed fines are for Amfg > and 
Am^3 < cases, respectively, and the dotted and the dash-dotted lines correspond to the 
same signs of Amf3 as above but with matter effect switched off. The mixing parameters 
are fixed as Amfg = ±3 x IQ-^ eV^, sin^ 2^23 = 1-0, Am?2 = 5 x IQ-^ eV^, sin^ 2^12 = 0.8, 
sin^20i3 = 0.05. We take pYg = 1.4 g/cm^ where p is the matter density and Yg is the 
electron fraction. 



As one might have suspected by looking in Fig.l the CP trajectory is elliptic. It 
is easy to prove it in the vacuum case and it can be also shown that it is the case 
in a good approximation for oscillations in matter. The proof of the statements will 
be given in Appendix. In particular, it is shown that the major (or minor) axis is 
always at 45 degree in the vacuum case. 

What does CP trajectory diagram actually represent? To answer the question, 
we start with the vacuum case and first concentrate on Fig.la-c. In these cases, the 
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lengths of the major and the minor axes are measures for the coefficients of sin 5 
and cos 5, respectively, in the oscillation probability. (The same statement holds for 
Fig. Id if the major and the minor axes are interchanged.) One can readily understand 
this statement by looking at Eq. (|8.4|) in Appendix. 

The two trajectories of positive and negative Am^g which are represented by 
dotted and the dash-dotted lines, respectively, are almost degenerate. One notices 
that the approximate degeneracy is between 6 and ti — 5 (mod. 2 tt) cases as we 
anticipated in discussions in the previous section. The two trajectories slightly split 
mainly due to the second term in Eq. (|2.3| ) and possibly by higher order terms of 



In matter, the CP trajectories of neutrinos and antineutrinos split; the former 
moves to downward-right (Am^g > 0) and the latter to upward-left (Am^g < 0). 
In fact, one can explicitly demonstrate that the matter effect is the cause of the 
departure of the two trajectories by artificially increasing the matter potential a. In 
Fig. 2, presented are the results of the same computations as in Fig. 1 but with 
a factor of 2 (artificially) larger matter effect. One observes that the neutrino and 
the antineutrino trajectories became more separated along the direction of major 
(minor) axis on bi-probability plane in Fig.2a-c (Fig.2d). Therefore, the degree of 
nonoverlapping of neutrino and antineutrino trajectories gives almost purely the 
measure for the matter effect. 

It is important to notice, by comparing the cases in vacuum and in matter in each 
figure, that lengths of the major and the minor axes, which measure the effect of sin 6 
term, the genuine CP violation, and the cos 6 term, respectively, barely change by the 
matter effect. It is nothing but the feature which we expect from the perturbative 
formula in Eqs. ( |2.3| ) and (P3|). 

Thus, the matter effect lifts the two-fold degeneracy of S and n — S cases dis- 
cussed in Sec. 2, as shown in Fig. 1. Notice, however, that there still remains two 
crossing points of the two trajectories corresponding to positive and negative Amfg 
as indicated in Fig. 1. It means that if we are unlucky so that the true value of 6 
is close to the crossing point, then we will still have two- fold degeneracy in determi- 
nation of 6. Since we expect only a modest statistics even with a huge detector in 
experiments for measuring 6, the region suffered by the ambiguity problem may not 
be so small, unless nature was so kind to tune her parameters to produce maximal 
CP violating effects. We will discuss in the next section how to resolve the remaining 
two-fold ambiguity. 

We emphasize that it is one of the nicest features of the CP trajectory diagram 
that the CP violating as well as conserving effects due to the Kobayashi-Maskawa 
phase, and the earth matter effect are pictorially displayed separately in a single 
diagram. 

The readers may worry about the apparently intricate features of the CP trajec- 
tory diagram indicated in Figs. 1 and 2. They include: 

(Fl) the feature that the "chirality" of the trajectory depends upon the sign of 
Am^g. Here, what we mean by "chirality" is how a trajectory winds, clockwise or 
counter-clockwise, as 6 varies from to 2tt. 




which are not taken into account in Eq. ( p.3|) . 
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Figure 2: The same as in Fig.0 but for pY^ = 2.8 g/cm^. 



(F2) the dependence on the neutrino energy; between energies of 0.5 GeV and 1 GeV 
the positions of 5 = and 6 = tt are exchanged completely. At the same time, the 
"chirality" of the trajectory also flips from 0.5 GeV to 1 GeV. 

The explanation of (Fl) and (F2) are in fact very simple. 

(Al) Since the 5-development of the trajectory is uniform, it suffices to discuss the 
behavior of the trajectory at around 6 = 0. When 6 increases from zero, the sin 6 term 
in the oscillation probability ( p.3| ) decreases (increases) in neutrino (antineutrino) 
channel. This means that the movement of the trajectory is toward upward-left 
direction for both Am^g > and Am^g < cases. 

The trajectory winds in opposite way for positive and negative Am^g because 
the point 5 = is located at near (far) side from the origin for Am^g > (Amfg < 0) 
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case, due to the cos 5 term in the probabihty ( p.3|) , at E = 0.5 GeV as in Fig. la. It 
means that the trajectory winds clockwise for Am^g > and winds counter-clockwise 
for Am^3 < 0. 

(A2) When neutrino energy increases from E = 0.5 GeV to 1 GeV, sin (— 
changes sign from negative to positive at L = 300 km. Then, the cos 6 term flips the 
sign and the 6 = point jumps to the far (near) side of the trajectory in Am^g > 
(Amjg < 0) case. Then, the CP trajectories winds to the opposite directions with 
those at E = 0.5 GeV. 



4. Principle of choosing beam energies for long-baseline neu- 
trino oscillation experiments 

The intriguing features of the CP trajectory diagram mentioned above implies a new 
principle for determining neutrino beam energies for given oscillation parameters and 
the baseline length. 

To avoid a confusion in notations that might occur when a major and a minor 
axes switch with each other depending upon the parameters we use in this section 
the terms the "radial thickness" and the "polar thickness" of contours. In Fig.la-c 
(Fig. Id), the polar thickness denotes the length of major (minor) axis, while the 
radial thickness implies that of minor (major) axis. 

As explained in the last section, the radial and the polar thickness of the CP 
trajectory diagram are the measure for cos 5 and sin 5 terms in the oscillation proba- 
bility, respectively. Now we try to maximize these thickness by tuning experimental 
parameters for a given set of mixing parameters. It is of course important to have 
a large sin 5 term because it is the signal for genuine CP violation. On the other 
hand, to maximize the radial thickness is to help in resolving the problem of two-fold 
ambiguity which was discussed in the last two sections. We believe that the latter 
is important in view of uncertainties due to statistical and systematic errors which 
would inevitably exist in any experiments. If we make a choice of, for example, 
E = 750 MeV for L = 300 km, then the contour shrinks to be approximately one- 
dimensional (see Fig. 6b) and there is no way to resolve the two-fold ambiguity p5 



We note that in Eq. (2^) the coefficients of the cos 5 and the sin 5 terms are 



proportional to xsinx and a;sin^(|), respectively, where x = — Therefore, a 
maximum of polar thickness, i.e., a maximum of absolute value of sin 5 term, occurs 
at half value of E/L of the corresponding maximum of radial thickness, coefficient 
of the cos S term. Namely, maximization of these terms implies: 



E \ f L \ ( Am 



Thus, maxima of polar thickness occurs at relatively low energies, (-E)sin(5 = 620 
MeV at L = 295 km and {E)sin5 = 1.5 GeV at L = 730 km, in agreement with 
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the conventional wisdom that CP violation effects are maximal at low energies. We, 
however, also emphasize that extremization of cos 6 term requires about twice larger 
values of beam energy for a fixed baseline distance. It is one of the reasons why we 
were led to the examination of a little bit of higher energies compared with that in 



Ref. |T6[, E = 0.5-2.0 GeV region, in this paper. The energies chosen in Fig.l turned 
out to be in the "right range" in compromising the requirements of maximizing the 
polar and the radial thicknesses of the CP trajectories. 

The readers might wonder the possibility that large errors due to statistical and 
systematic uncertainties completely invalidate our principle of optimization of beam 
energies. We will demonstrate in Sec. 6 that it does not occur at least for certain 
range of reasonable oscillation and experimental parameters. 



5. Resolving two- fold ambiguity in determination of 6 

We have shown in the previous sections that there exists a two-fold ambiguity in 
determination of CP violating phase 6 due to our ignorance of the sign of Am^g. We 
discuss in this section the problem of how to resolve the ambiguity. 

5.1 Chance for simultaneous measurement of S and the sign of Am^g in situ 

Our foregoing analyses of the two-fold ambiguity have revealed the intriguing possi- 
bility that, if we are lucky, an in situ simultaneous determination of S and the sign 
of Amfg may be possible even at relatively short baseline as L = 300 km. Namely, 
the earth matter effect is comfortably large to split the two trajectories correspond- 
ing to positive and negative Am^g such that the both quantities can be measured 
simultaneously in a certain range of 6. Namely, if the angle 6 is in the third or the 
fourth quadrants for positive Amfg (the normal mass hierarchy), or if S is in the first 
or the second quadrants for negative Am^g (the inverted mass hierarchy), then mea- 
surement of Pii'fj, Ve) and P{i'^ — ^ v^,) in neutrino and antineutrino experiments 
can determine both quantities simultaneously. 

The statement just made above is a conservative one and the range of lucky 
determination of 5 without ambiguity may extend to wider region, as one notices in 
Fig. 1. However, how wide is the region depends upon the mixing parameters as 
well as experimental uncertainties. We describe an attempt toward quantifying it in 
Sec. 6 by computing numbers of events. 

Our above observation sharply contrasts with the conventional belief that very 
long baseline as > 1, 000 km is required for the determination of the sign of Am^g, 
and open the door to a simultaneous measurement of 5 and the sign of Am^g in situ 
at relatively short (~ 300 km or so) baseline. 

We will further pursue the possibility of simultaneous measurement of 5 and the 
sign of Amfg in a single experiment which is valid for full range of 5 by considering 
longer baseline in Sec. 6. 

5.2 Two-detector method 

However, nature may not be so kind. Namely, if the true value of 5 is within the 
experimental uncertainties to one of the two crossing points of the two trajectories. 
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Figure 3: Asymmetry of the probability ratio defined in Eq. (5.2) in the text for (a) 5 = 
and (b) 5 = tt/2. The mixing parameters are chosen as the same with those of Fig. 1. 



we have a two-fold ambiguity in the determination of 5. In the worst case in which 
6 is really close to the crossing points, one cannot resolve the two-fold ambiguity no 
matter how accurate were the measurement. 

Then, the question is how to resolve the two-fold ambiguity. In this subsection, 
we discuss the two-detector method in order to resolve the two-fold ambiguity in 
determination of the CP violating angle 6, in case it remains in single-detector ex- 
periments. The multiple detector method has been proposed by various authors for 



different physical motivations ||34 |. 



We first illuminate that the two-detector option is naturally motivated by the 
nature of the phenomenon of neutrino oscillation itself. Let us first define the ratio of 
the suitably energy averaged appearance probabilities for neutrino and anti-neutrino. 
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R{P), as follows, 



We note that if there is no CP violating phase and matter effect, R{P) does neither 
depend on the baseline nor on the average neutrino energy as long as energy distri- 
butions of neutrino and anti-neutrino are the same. Some dependence of R{P) on 
the baseline can indicate either genuine CP violating effect or matter effect, or both. 
However, by simply looking at R{P), it may be difficult to separate these two effects. 

In order to see for which energy and baseline the matter effect could be important, 
independent of the magnitude of CP phase, let us define the asymmetry of the ratio 
R{P) as follows, 

,(p^-.., RjP-^ ^^Is > 0) - RjP; A^i3 < 0) 

A large value of the asymmetry implies that the matter effect is enhanced relative to 
the vacuum effect. In fact, if there is no matter effect, from the expression of vacuum 
probability ( p.3| ), this quantity is expected to be small for any values of 6. 

In Fig. 3 plotted is the asymmetry A{R) for two typical values of CP violating 
angle, 6 = and 5 = |. The mixing parameters are chosen as the same with those 
of Fig. 1. There are several characteristic features in the figures. First of all, as 
expected, the asymmetry depends very weakly on 6. The asymmetry tends to be 
positive, indicating that Am^g > case is enhanced compared to Amfg < case. 
There are two distinct regions in which the asymmetry A{R) is large and positive. 
At E ~ 1 GeV, they are at L ~ 600 - 700 km, and at 1000 - 1500 km. 

We have repeated the similar computation to obtain the CP trajectory diagram 
for the path lengths L = 700 and 1000 km, and energies E = 0.5, 1.0, 1.5, and 2 
GeV. The results are presented in Fig. 4 and 5. By looking into Figs. 1, 4 and 5, we 
immediately observe several distinct features: 

(1) There are cases with longer baselines in which the neutrino and the antineutrino 
trajectory do not intersect; E = 1.0, 1.5, 2.0 GeV at L = 700 km, and E = 1.5, 2.0 
GeV at L = 1000 km. We note that this is consistent with what we presented in 
Fig. 3 where we can see that these values of experimental parameter fall in the 
region of strong asymmetry defined in Eq. ( ^.2|) . It resolves the two-fold degen- 
eracy we talked about. Of course, by focusing these experimental parameters, one 
can determine 6 and the sign of Am^g simultaneously, provided that statistical and 
systematic uncertainties are small enough. It raises an important possibility that 
simultaneous determination of 6 and the sign of Am^g can be done in full range 
of 6 with optimal distance of L ~ 700 km. We will further discuss the possibility 
including experimental uncertainties in the next section. 

(2) The path-length dependence of the trajectory diagram is not always smooth. For 
example, the feature displayed in E = 1 GeV at L = 1000 km case does not fall into 
a smooth extrapolation of the behavior oi L = 295 and 700 km at the same energy. 
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Figure 4: The same as in Fig.0 but for L = 700 km. 



(3) Some curious behavior is observed in the cases, E = 0.5 GeV at L = 700 
and 1000 km, where the matter probabihties have similar behavior with those of 
vacuum oscillation, which seems to be also consistent with the results in Fig. 3. 
Most probably, it is a new phenomena, not a remnant of the vacuum mimicking 
phenomenon, whose interpretation is not known at the moment. 



A tentative conclusion before examining the numbers of events is that by using 
the second detector placed at distances of 700 km or 1000 km, the two-fold ambiguity 
will be resolved if it remains in shorter baseline JHF-type experiments. 
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6. Long-baseline neutrino oscillation experiments with super- 
beams; some concrete examples 

In discussions in the foregoing sections we have concentrated on illuminating global 
structure of the neutrino oscillations in the region relevant for low energy experi- 
ments, and did not pay enough attention to the statistical and the systematic un- 
certainties. Of course, they are of key importance in judging what is the optimal 
design of the experiments. On the other hand, it is not quite possible to determine 
at this moment which is the optimal design because we do not know the value of 
key parameters, most importantly 9i^. Therefore, our discussion in this paper is 
inevitably restricted to the one that may be called, at best, case studies. 

In this section we first try to estimate the numbers of events including signal and 
the background by taking the same typical mixing parameters as the one we used 
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before, which correspond to the LMA solution of the solar neutrino problem. We do 
this by taking the concrete examples of upgraded long-baseline neutrino oscillation 
experiment which utilizes a conventional neutrino superbeam. For definiteness, we 
assume a 4 MW proton beam which is planned to be constructed in the Phase II 
of the JHF project, and consider 1 Mton Super-Kamiokande type water Cherenkov 
detector whose fiducial volume is assumed to be 0.9 Mton. Such detector is already 
discussed by the Super-Kamiokande group who mainly motivated by an extensive 
search for proton decay under the name of Hyper-Kamiokande project |3^. The goal 
of our study in the present section is not to develop the experimental proposal which 
is ready to submit, but to obtain a feeling on what are the promising possibilities 
which deserve further detailed studies, possibly by joint collaboration by theorists 
and experimentalists. 

We discuss three options, 

(i) Single detector at L = 300 km. 

(ii) Single detector at L = 700/1000 km. 

(iii) Two detectors at L = 300 and 700/1000 km. 

In this work, we will use two types of different neutrino beams with quasi- 
monochromatic energy spectrum calculated by the JHF-SK Neutrino Working Group 
p9[] , which are (i) narrow band (NB) beam and (ii) off axis (OA) beam. NB beam is 
made by pions with particular choice of momentum selected by dipole magnet placed 
between two horns in the wide band beam configuration. We use in our analysis two 
different options of NB beam characterized by pion momenta of 2 and 3 GeV, which 
peaks at ~ 1 GeV and 1.4 GeV, respectively. OA beam is an another way of making 
quasi- monochromatic neutrino beam proposed in Beavis et al. in [Q. It is obtained 
by slightly (a few degrees) displaced the direction of axis of wide band beam from 
the far detector direction. For this type of beam, we use the ones obtained with 
off axis of 3 and 2 degrees, which peaks at ~ 0.5 GeV and 0.8 GeV, respectively. 
While we adopt as the basic parameters the ones calculated for the JHF neutrino 
experiment, we do hope that the results of our calculation are illuminative enough 
for future projects on other continents as well. 



6.1 Method for calculation of number of events 

Before proceeding to physics discussions we have to explain first how we calculate 
expected numbers of events. In our computation, we take into account both signal 
and background, which are calculated by the way as explained below. 

Signal consists of the contributions from charged current (CC) Vf, (coming from 
oscillation) interactions, which are classified into quasi-elastic (qe); i/ + A^ — 
i + A^', one pion production (Ivr); z/ + A^^£ + A^' + 7r, multi pion production (mvr); 
z/ + A^— i>£ + A^' + nvr, and coherent pion (cvr) production; v Q ^ v O + 7r+, 
reactions. In our computation, we define the signal coming from the CC interactions 
in the expected number of events, N^ig, as follows, 

AT,,, ^ T dEa'cc{E)<Pl^{E)P,,{E)t{E), (6.1) 

i=ge,l7r,m7r,C7r 
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where and a'^cc{E) are the the number of target and the CC cross section, 
respectively, for i-th reaction process, T is the exposure time, (pl^^E) is i/^ neutrino 
flux at the detector site in the absence of oscillation as a function of neutrino energy, 
P^e{E) is the z/^ — > Ue oscillation probability, and e{E) is the detection efficiency for 
the e-like events, which was taken from Ref. |]30|. See Refs. [0, for the neutrino 
cross sections we used in this work. 

We also take into account possible background which are coming from 7r° pro- 
duced in the NC and the CC interactions, e//i misidentification, and u^. contamination 
in the original z/„ beam, where the dominant ones come from 7r° produced in NC re- 



actions as well as contamination. Following Ref. |^ , we express background event 
as 

Nbg = N^S + ^bg + ^Ig™, (6-2) 
where N^q, N^q and N^^q^ are backgrounds coming from the CC, the NC inter- 
actions and the z/g contamination, respectively. We compute N^q in a similar way 
as Nsig in Eq. ( |6.1D but P^e{E) replaced by P^^{E) and e{E) replaced by a con- 
stant reduction efficiency ~ 0.017 %, estimated from Table 2 in Ref. [0, as an 
approximation. 

Similar to the CC interactions, contributions from NC ones include the following 
reactions: elastic (es); v+N — > v+N , one pion production; v + N u + N' + t[, multi 
pion production; v + N-^v + N' + mr, and coherent pion production; u O 
V O + 7r°. We can express as 

N^E^riT E Ni.JdEa^r,c{E)<Pl{E), (6.3) 

i=es ,l7v ,rmT ,ciT 

where rj is the reduction efficiency, approximated as constant, which is about 0.22 % 



of the total NC events, estimated from Table 2 in Ref. |^ . 

The contribution from the beam z/g contamination can be estimated in a similar 
way as in eqs. ( |6.1D and ( |6.3| ) but [E) replaced by 0°^ {E) , the z/g flux which exists 
in the z/^ flux in the absence of oscillation and P^e{E) by Pee{E) in Eq. ( p.l|) . 

We deflne the expected number of "e-like" events in —>■ z/g channel as the sum 
of signal and background as, 

iV(e-) = N,,g + Nbg. (6.4) 

For anti- neutrino channel, z/^ —>■ Ue, we compute the expected number of e+-like 
events in the same way properly replacing neutrino flux as well as cross sections by 
that of anti-neutrino. 

In Table 1 we present expected number of events without (which means to set 
e = rj =100 % in Eqs.( |6.1| ) and ( |6.3D ) and with detection efficiency for JHF neutrino 



OA 2 degree beam with baseline 295 km and 0.77 MW beam power assuming Super- 
Kamiokande detector (22.5 kton), corresponding to the first phase configuration of 
the JHF neutrino project |31[|, with exposure of 5 years. For comparison we also 
present the numbers found in new JHF LOI The numbers we obtained are 



rather similar to those computed by the JHF working group, allowing us to have 
confidence on that our computations are accurate enough for our purpose. 
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CC 




Beam Vf. 


Oscillated Vf. 


This work (before reduction) 


10707 


3840 


270 


297 


This work (after reduction) 


1.8 


8.8 


10.3 


142.7 


JHF Phase I (before reduction) 


10714 


4080 


292 


302 


JHF Phase I (after reduction) 


1.8 


9.3 


11.1 


123.2 



Table 1: Comparison of our expected number of events with that of the JHF working 
group for OA 2 degree beam option [^] with basehne 295 km, 0.77 MW beam power 
and the exposure of 5 years assuming the SK detector (22.5 kton). For the calculation 
of oscillated fg signal, Amfg = 3 x 10~^ eV^, sin^ 26*23 = 1-0, sin^ 26*13 = 0.1, and other 
mixing parameters are set to be zero (Amf2 — > 0, sin^ 2^12 — > 0) and matter effect was 
neglected. For each case, upper and lower number indicate the expected number of events 
before reduction (or with 100 % detection efficiency) and after reduction (with reduction 
and/or detection efficiency). 



6.2 Single detector at L = 300 km 

This is essentially identical with the upgraded JHF neutrino experiment in its Phase 
IL While its detailed study is underway in the working group, we try to make some 
suggestions here which might add (we hope) their final proposal some useful ingre- 
dients. 

In order to have some feelings about the expected number of events for various 
JHF neutrino beam options, we show in Table 2 the numbers of electron appearance 
events assuming 100 % Vfj, ^ Ve conversion, for 1 megaton SK type detector (with 
fiducial volume of 0.9 Mton) with baseline L = 295 km, and 4 MW beam power as 
planned in the Phase H of the JHF neutrino experiment with an exposure of 1 year. 
We also show the background which is essentially constant and does not depend 
on neutrino oscillation. One can estimate the expected number of events in the 
presence of neutrino oscillation from the numbers in the table by simply multiplying 
oscillation probability which is properly averaged over the cross section, neutrino 
energy spectrum, detection efficiency, etc. 



As we emphasized at the beginning of this section the statistics heavily depends 
on the value of ^13, on which we have no clue at present apart from the CHOOZ 



constraint |T0[- Hence, it is difficult to judge if the experiment that we are discussing 
is feasible or not at this moment. Therefore, we restrict ourselves in this paper to a 
particular value of ^13, sin^ 2^13 = 0.05, which is about a half of the CHOOZ bound. 

Now we present in Fig. 6 the CP trajectory contour written on the plane spanned 
by numbers of events (not the probabilities), A^(e~) and N{e~^), including the back- 
ground assuming 2 and 6 years running for neutrino and anti-neutrino channels, 
respectively. The Fig. 6a is for OA 3 degree beam which peaks at £^ ~ 0.5 GeV, 
while Fig. 6b is for OA 2 degree beam which peaks at ~ 0.8 GeV. Fig. 6c and 
6d are for NB 2 and 3 GeV beams, whose energy spectra peaks at ~ 1 GeV and 
1.4 GeV, respectively. We note that only for NB 3 GeV beam, we assume, as an 
approximation, anti-neutrino flux is the same as that of neutrino since the flux for 
anti- neutrino was not available for this particular beam option [^]. Dotted circles 
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Narrow Band Beam 


Off Axis Beam 




2 GeV 


3 GeV 


3 deg. 


2 deg. 


Peak Energy (GeV) 


~ 1.0 


~ 1.4 


~ 0.55 


~ 0.75 


V (signal) 


1.8 X 10^ 


2.7 X 10^ 


1.6 X 10^ 


4.3 X 10^ 


V (background) 


1.7 X 10^ 


2.5 X 10^ 


3.5 X 10^ 


8.2 X 10^ 


V (signal) 


6.3 X 10^ 


1.2 X 10^ 


4.6 X 10^ 


1.4 X 10^ 


V (background) 


0.6 X 10^ 


1.2 X 10^ 


1.1 X 10^ 


2.7 X 10^ 



Table 2: Expected numbers of e~ and e^-like events in neutrino and anti- neutrino 
channels for some possible JHF neutrino beam options we use in our work, assuming 100 
% conversion of — fg and z/^ — > i/e, respectively. We take a 1 Megaton SK-type water 
Cherenkov detector (assuming fiducial volume of 0.9 Mton) with baseline distance of 295 
km, 4 MW beam power (corresponding to the Phase II of JHF neutrino project), and 
exposure of 1 year. Signal as well as background events are presented separately in the 
table. 

in Fig. 6 indicate uncertainty which corresponds to 3 cr, where only the statistical 
error is taken into account. Unfortunately, we have no way of estimating expected 
systematic errors in the super- JHF experiments. 

As we noticed in Sec. 3 there is an inversion phenomenon in the diagram. 
Namely, for a given shape of the contour, the low probability branch (cos 5 > 
region) in Fig. 6a is mapped into the high probability branch (cos 5 < region) in 
Fig. 6b in Am^g > case, and vice versa in Am^g < case. 

We point out that the feature can be used as a method of identifying the value of 
5, if otherwise ambiguous in its measurement. Suppose that they run first at = 1 
GeV and obtained the result which tends to prefer the high probability branch, but 
they were not confident because the high and the low probability branches are not 
so well separated. While they can just continue to run to increase statistics with the 
same energy in such circumstance, an alternative way (and better way, we believe) 
is to run at lower energies, ~ 0.5 GcV. If the parameter is really in the high (low) 
probability branch sX E 1.0 GeV, then it must jump down (up) to the low (high) 
probability branch at £■ ~ 0.5 GeV. 

Let us consider, for example, the case of positive Am^g. If 5 ~ | (high probability 
branch), the jump is from P{y) ~ 2.05 % at 1 GeV to P{y) ~ 1.3 % at 0.5 GeV, 
as one can see in Fig. 1. On the other hand, if 5 ~ ^ (low probability branch) 
the jump is from P{v) ~ 1.85% at 1 GcV to P{v) ~ 1.65 % at 0.5 GeV. Thus, a 
greater downward jump of about 37 % is predicted in the case of 5 ~ ^, as compared 
to a modest ~ 10 % decrease in the 5 ~ ^ case. Thus, running at high and low 
energies would greatly help to distinguish the high and the low probability branches, 
or cos (5 > and cos 5 < cases, if the energies are tuned to have a large radial 
thickness. Unfortunately, wc cannot resolve the two-fold ambiguity by using this 
technique, because the high probability branch jumps to low probability branch, or 
vise veasa, both in Amfg > and Amfg < cases simultaneously. 

It should be emphasized that, if we are lucky, we will be able to determine the 
CP violating angle 5 and the sign of Amfg simultaneously. It is the case with highest 



19 



4000 



(a) OA 3deg. 



Am3(>0 
Am^3i<0 



3000 







2000 




1000 ^ 
2000 

3200 



2800 



2400 



2000 



1600 



1200 



3000 



4000 



9000 



(b) OA 2deg. 





8000 
7000 - 
6000 - 

5000 : 

4000 r 

5000 5000 6000 7000 8000 9000 10000 

2800 I 

(d) NB 3GeV 

2600 - 



2400 
2200 
2000 
1800 
1600 




1800 2200 2600 3000 3400 3800 1600 1800 2000 2200 2400 2600 2800 

N(e") N(e") 

Figure 6: CP trajectory in the N{e~) — N(e~^) plane for off axis and narrow band beam 
for L = 295 km. Dotted circles indicate 3 a errors (only statistical ones) . 



statistics among the cases considered in this paper, and hence it is hkely to have 
resolving power of the both quantities. By luckiness, we mean very roughly that if 
the angle 6 is in the third or the fourth quadrants in the case of normal mass hierarchy 
(Amjg > 0), and is the first or the second quadrants in the case of inverted mass 
hierarchy (Am^g < 0), respectively. In this case, the experiment will reveal by high 
statistics data taking the CP violating angle and the sign of Am^g simultaneously. 
The accuracy of the measurement of 6 would be somewhere in 20-30 degrees, if the 
systematic errors are not too sizable. 

6.3 Single detector at L = 700/1000 km 

We next consider the option of single detector at L = 700/1000 km. As we saw in Sec. 
3 the detector at L = 700/1000 km shows a qualitatively different characteristics; 
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Figure 7: CP trajectory in the N{e~) — N(e~^) plane for off axis and narrow band beam 
for L = 500 km or 700 km. Dotted circles indicate 3 a errors (only statistical ones). Mixing 
parameters are fixed to be the same as in Fig. 1. 




the neutrino and the antineutrino trajectory becomes nonover lapping for neutrino 
energies of 1 — 2 GeV, aUowing simultaneous determination of S and the sign of Am^g 
for any values of 6, if the statistics is high enough. 

In Fig. 7 we present the CP trajectory contours written on the N{e~) — N{e~^) 
plane with path length of (a, b) 500 km and (c, d) 700 km. The former is for 
comparison, in particular to represent the correlation between the matter effect and 
the number of events as a function of baseline length. Since there is no version of 
off-axis beam with peak energies of > 1 GeV available to us, we use NB 3 GeV beam 
which peaks at around E — lA GeV. It appears that off-axis beams are more intense 
than narrow-band beams with similar peak energies by a factor of 2-3. If we are able 
to design OA beam with peak energy of > 1.4 GeV, we would gain a factor of 3 in 
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number of events. If it is the case, we may scale the abscissa and the ordinate by 
factor of 3 and the radii of the error circle are reduced by factor of y/S. Of course, to 
have off-axis beam with peak energies of > 1 GeV one should probably start with 
a completely different proton beam design with higher energies, whose discussion is 
far beyond the scope of this paper. 

We note that for NB 3 GeV the number of events is ~ 500 for 2 (6) years of 
running in neutrino (antineutrino) channel at L = 700 km, and ~ 2000 for L = 500 
km. They are certainly not enormous but are not too small either. If the beam 
normalization is properly known it would be possible to discover the existence of 
leptonic CP violation and at the same time to determine the sign of Am^g. 

Clearly, the conclusion depends on many unknown factors, the values of mixing 
parameters, ^13, Arn^g^ ^ ^Iso one experimental conditions, uncertainties in 
absolute flux normalization, cross sections, and efficiency background rejection. 



6.4 Two detectors at L = 300 and 700/1000 km 

Suppose that we have a bad luck in the experiment with L = 300 km by having 5 
in the alternative regions from that we have mentioned in Subsec. 6.2. Namely, if 
sin 5 > in the case of normal mass hierarchy and sin 5 < in the case of inverted 
mass hierarchy, we will be able to determine (assuming enough statistics) 6 but only 
with modulo two-fold ambiguity. While it is already a great achievement, it may be 
better if we have a way to resolve the ambiguity. Motivated by the consideration 
in the preceding section, it is natural to consider the two-detector option, one at 
L = 300 km and the other at L = 700/1000 km. While it is indeed an "expensive 
option" which uses the two megaton detectors (or possibly one 1 Mton and one ~ 
100 kton iron calorimeter) it might not be an unrealistic one in view of the proposal 
of detectors either in Korea |3^, or in Beijin . 

In this two-detector option, the requirement of statistics can be relaxed because 
the second detector is effectively only for the determination of the sign of Am^g. In 
view of Fig. 7 a narrow band beam with energy a bit higher than NB 3 GeV beam 
may be able to do the job, without multiplying 3 in flux normalization 



7. CP trajectory diagram for neutrino factories 

Finally, we briefly address in this section the features of CP trajectory diagram for 
physical parameters which are appropriate for neutrino factories. It is to illustrate 
enormous difference in the features of the CP phase-matter interplay between the 
situations of neutrino factory and the low energy superbeam. We neither intend to 
make full comparison of these two different strategies nor try to argue which is the 
better way to measure CP violation. 

In Fig. 8, presented is the CP trajectory diagram for L = 3000 km for varying 
neutrino energies from 10 GeV to 50 GeV. Only the trajectories with matter effect 
are plotted for both sign of Am^g. The matter density is taken as pYf, = 2 g/cm^ 
corresponding to upper mantle region. Unlike diagrams for low energy beams pre- 
sented in the preceding sections the two trajectories corresponding to positive and 
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Figure 8: The CP trajectory diagram in bi-probability plane for L = 3000 km and much 
higher neutrino energies E = 10—50 GeV which correspond to so called "Neutrino Factory" 
situation. The mixing parameters are fixed to be the same as in Fig. 1 except that we take 
pYe = 2.0 g/cm^. 



negative Amfg are far apart with each other. It leads to the well known fact that the 
matter effect dominates over the CP violating effect due to the Kobayashi-Maskawa 
phase. 

The good news is that in the neutrino factory situation there is of course no 
such ambiguity as the one we discussed in this paper. (However, there arises an 
another ambiguity in this case, if the value of Oi^ is not known, as pointed out by 
Burguet-Castell et al. |28|, as mentioned in |3||.) On the other hand, the problem of 



matter effect contamination is severer. It appears to us that we need further studies 
to clearly separate genuine CP violating effect from the matter effect in neutrino 
factories. 
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8. Conclusion 



We have discussed the features of the neutrino oscillations which are relevant for 
the experiments with conventional low energy (typically ~ 1 GeV) superbeams with 
mediumly-long baseline ( < 1000 km) which aim at measuring leptonic CP viola- 
tion. We have assumed, for the purpose of estimating number of events and statis- 
tical uncertainties, a supermassive 1 megaton water Cherenkov detector of Super- 
Kamiokande type, and 4 MW beam power as planned in the Phase II of the JHF 
neutrino experiment. 

We started with focusing on the problem of interplay between genuine CP vio- 
lation and the matter effect. While it is a widely discussed topics in the literature, 
we have uncovered an interesting new feature; the matter effect helps. Namely, the 
matter effect in such mediumly-long baseline experiments is comfortably large so 
that it helps to resolve the (approximate) two-fold ambiguity which would exist, in 
the absense of the effect, in the vacuum oscillation probability. 

To elucidate the features of the CP phase-matter interplay, we have introduced 
a new powerful tool called the "CP trajectory diagram in bi-probability space" . The 
diagram allows us to represent pictorially the following three effects separately in 
a single diagram in a compact form: 1) the genuine CP violating effect due to the 
leptonic Kobayashi-Maskawa phase 5, or sin S, the CP odd term, 2) the effect due to 
the CP even term, cos 5, and 3) the fake CP violating effect due to the earth matter. 
As discussed in the text, the effect of the CP odd term is characterized by the "polar 
thickness" whereas the effect of the CP even term is characterized by the "radial 
thickness" of the trajectory contour (see the beginning of Sec. 4 for the definitions 
of "polar" and "radial" thinkness). On the other hand, the earth matter effect is 
characterized by the distance of separation between the two trajectory contours with 
different sign of Aml^. 

By using such diagram, we have observed that there is a two-fold ambiguity in 
the determination of the CP violating angle 5, if we do not know a priori the sign 
of Am^g. The ambiguity is shown to be a remnant of the approximate symmetry 
possessed by the vacuum oscillation probability (as mentioned above) under simul- 
taneous transformation of 5 — >• tt — 5 and Amf^ — >• — Amfg. We have discussed the 
principle of "maximal fatness" of the CP trajectory diagram to find optimal exper- 
imental parameters in order to have a large CP violation and at the same time to 
resolve such two-fold ambiguity. 

By these considerations we are naturally invited to consider the possibility of 
measuring simultaneously the CP phase 6 as well as the sign of Am^g. We have 
discovered the enlighting possibility that such simultaneous measurement can be done 
with relatively short baseline, L ~ 300 km assuming 1 megaton water Cherenkov 
detector, as in the phase 11 of JHF-Kamioka neutrino experiment. It works, however, 
only under the condition of nature's kind setting of the parameter into the region 
sin 6 ■ Amjg < 0. 

We also described a way of resolving the degeneracy in the case of unresolved high 
and the low probability branches, or cos 5 > and cos 5 < region for Am^g > 
case, and vice versa in Am^g < case. By running at high (~ 1 GeV) and low 
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(~ 0.5 GeV) beam energies the high and low probabihty branches are interchanged 
with each other, and this effect may help in identifying 6 while it cannot resolve the 
two-fold ambiguity. 

Furthermore, we have uncovered an interesting possibility that simultaneous 
measurement of 6 and the sign of Am^g can be done in situ for the whole region 
of 5 (0 < (5 < 27r) in the experiments with longer baseline distance, L ~ 700 km 
(corresponding to the CERN to Gran Sasso and/or Fermilab to Soudan-2 distances), 
again with a megaton water Cherenkov detector. 

We also briefly discussed how the CP trajectory diagram look like for the neutrino 
factory situation and observed that the two CP trajectory with different signs of 
Am^3 are well separated owing to the fact that the matter effect is dominant. It 
is in sharp contrast with the case of low energy superbeam, the main subject of 
this paper. A clear way of separating the CP phase effect from the matter effect is 
naturally called for. 

We hope that our discussions in this paper would help in designing the future 
neutrino oscillation experiments. 
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APPENDIX 

We give a proof that the CP trajectory is elliptic in this Appendix. It can be 
shown under the assumption of mass (Am^) hierarchy and the adiabatic approxima- 
tion that the neutrino and the antineutrino oscillation probabilities can be written 
in the following forms: 



P{u) = Acos6 + Bsm6 + C 
Piv) = A cos 5 - 5 sin 5 + C* 



U) 
1.2) 



It is true in all the perturbative formula so far derived in Refs. |20, 21, 22, 23]. 
However, we should note that it is not true in the exact formula derived in . 
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Therefore, it is not an exact formula apart from the case in vacuum, but is a very 
good one in the case of hierarchical mass difference we are interested in. 

Once the above general form holds for a fixed arbitrary energy, it hold even if one 
take average over an arbitrary neutrino energy spectrum. Therefore, we implicitly 
imply P's and the coefficients in (|8.2| ) as averaged over an energy spectrum, though 
it is not indicated explicitly. Hence, the general form applies to the CP trajectory 
diagram on bi-number-of-event plane, with which we have dealt in Sec. 6. (In fact, 
this general form applies even if we take average over the neutrino path length.) 

By eliminating 6 we obtain the equation obeyed by -P(z^) and -P(z^) as 




In vacuum, the expression simplifies because A = A etc.: 

1, (8.4) 



\ ) \ 2B ) 

which implies that the minor (if A > B) oi major (if A < B) axes are always at 45 
degree. 

This completes the proof that the CP trajectory is elliptic. 
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